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Abstract 

Recent analytical and numerical work have shown that the spectrum of 
the random non-hermitean Hamiltonian on a ring which models the physics of 
vortex line pinning in superconductors is one dimensional. In the maximally non- 
hermitean limit, we give a simple "one-line" proof of this feature. We then study 
the spectral curves for various distributions of the random site energies. We find 
that a critical transition occurs when the average of the logarithm of the random 
site energy squared vanishes. For a large class of probability distributions of the 
site energies, we find that as the randomness increases the energy at which 
the localization-delocalization transition occurs increases, reaches a maximum, 
and then decreases. The Cauchy distribution studied previously in the litera- 
ture does not have this generic behavior. We determine 7ci, the critical value 
of the randomness at which "wings" first appear in the energy spectrum. For 
distributions, such as Cauchy, with infinitely long tails, we show that 7ci = 0+. 
We determine the density of eigenvalues on the wings for any probability distri- 
bution. We show that the localization length on the wings diverge generically 
as L{E) ~ ^E-Et\ ^® ^ approaches E^,. These results are all obtained in the 
maximally non-hermitean limit but for a generic class of probability distributions 
of the random site energies. 



1 Introduction 



Recently, there has been a great deal of work on non-hermitean matrices , operators 
0,0, and Hamiltonians To make this paper as self-contained as possible, we will 
begin with a brief review of the material relevant for our discussion. Non-hermitean 
operators appear in a broad class of problems in statistical physics describing diffusion 
and growth exemplified by the equation 

dP 

with Pj representing for example the concentration on site j or the population of 
species j. The w's denote various transition probabilities and growth rates. Writing 
©as 

(2, 

we see that the operator L is in general non-hermitean. 

A non-hermitean Hamiltonian inspired by the problem of vortex line pinning in 
superconductors has attracted the attention of a number of authors |^, ||, H, |Tl 
In one version, we are to study the eigenvalue problem 



with the periodic identification i + = i of site indices. Here the real numbers 
Wi are drawn independently from some probability distribution P iw) (which we will 
henceforth take to be even for simplicity.) This Hamiltonian describes a particle 
hopping on a ring, with its clockwise hopping amplitude different from its counter- 
clockwise hopping amplitude. On each site there is a random potential. The number 
of sites is understood to be tending to infinity. We will also take t and h to be 
positive for definiteness. (The hopping amplitude t can be scaled to 2 for instance but 
we will keep it for later convenience.) Note that the Hamiltonian is non-hermitean 
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for h non-zero and thus has complex eigenvalues. It is represented by a real non- 
symmetric matrix, with the reality implying that if E is an eigenvalue, then E* is 
also an eigenvalue. For P (w) even, for each particular realization {wi} of the random 
site energies, the realization {— w,} is equally likely to occur, and thus the spectrum 
of H averaged with P (w) has the additional symmetry E ^ — E. This eigenvalue 
problem is clearly a special case of (|lD. 

Without non-hermiticity {h = 0), all eigenvalues are of course real, and Anderson 
and collaborators ||T^ showed that all states are localized. Without impurities {wi = 
0), Bloch told us that the Hamiltonian is immediately solvable with the eigenvalues 

Er, = t cos - ih) , (n = 0,l,---,iV-l), (4) 

tracing out an ellipse. The corresponding wave functions ipj^^ ~ exp 2Tiinj /N are 
obviously extended. We are to study what happens in the presence of both non- 
hermiticity and the impurities. 




Figure 1: The eigenvalues of one particular realization of the Hamiltonian in (3) with 
= 300, t = 2, h = 0.5 and with the site energies drawn from the box distribution 
(13) with V = 1.8. 

With impurities present, two "wings," a right wing and a left wing, emerge out of 
the two ends of the ellipse. (See figure 1.) Some eigenvalues are now real. Evidently, 
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the "forks" where the two wings emerge out of the "elhpse-hke" curve represent a non- 
perturbative effect, and cannot be obtained by treating the impurities perturbatively. 
For P {w) even, we can focus on the right wing in hght of the symmetry mentioned 
above. 

Le Doussal |T3| and Hatano and Nelson emphasized that H has a special prop- 
erty, namely that by a (non- unitary) gauge transformation ^ ipj = e~^^ ipj, all 
the non-hermit icity can be concentrated on one arbitrarily chosen bond, on which the 
hopping amplitudes are | without changing the spectrum. Indeed, performing 

this transformation on @ we obtain the eigenvalue problem H ip = Eip with 







Hij = 




Hni = 




2 


HiN = 


^ p-Nh 


2 



(5) 

The non-hermit icity has been moved to be on only one particular link, which in (^) 
is between sites 1 and A^, but it is clear that by a suitable transformation we can put 
the hopping amplitude | on any link. 

This leads to an extraordinarily simple argument 0] that the states corresponding 
to complex eigenvalues are extended, that is, delocalized. Assume that an eigenstate 
oiif) oiH is localized around some site j. We can always gauge the non-hermiticity to a 
link which is located as far away from the site j as possible, where is exponentially 
small. Thus, if we cut the ring open at that link, the effect on the Schrodinger equation 
would be exponentially small, if the localization length of the state is smaller than 
h~^. The state does not know that the Hamiltonian is non-hermitean and so its 
eigenvalue must be real. Transforming back to if, we have a localized state whose 
eigenvalue is real. On the other hand, if h is larger than the inverse localization length 
of this state, we no longer have a normalized eigenstate of H . Hatano and Nelson 

argued that such states of H and extended states of H correspond to states of H 
with complex eigenvalues. They thus concluded that the states of H corresponding 



to complex eigenvalues are extended, that is, delocalized. 

Remarkably, non-hermitean localization theory is simpler in this respect than the 
standard hermitean localization theory of Anderson and others |jT^. To understand 
the localization transition, we have to study only the energy spectrum of H . 

One way to think about this problem is to imagine starting with Hq (that is, H 
with {wi\ set to zero) and then increasing the randomness (measured by a parameter 
7, say) and ask if there is a critical strength of randomness 7ci at which localized states 
first appear. There may also be another critical strength 7c at which all the extended 
states become localized. Obviously another way of thinking about the problem is to 
start with the hermitean random problem {h = 0) in which all states are localized. 
We then increase the non-hermiticity and ask about the critical strength hf. of the 
non-hermiticity at which extended states first appear. Similarly, we can ask, for given 
h and 7, the critical hopping amplitude tc needed for the states to become delocalized. 

As we have mentioned, a number of authors have considered (^ . In this paper, we 
will refer extensively to |jlO[ and [|lT] which we will denote by FZ and BZ respectively. 
Here, as in FZ and BZ, the emphasis is to understand analytically the essential physics 
involved. 

We will now describe some of the results from FZ and BZ, partly to fix the 
framework for our subsequent discussion. As is standard, we are to study the Green's 
function 

= (6) 

The averaged density of eigenvalues, defined by 

p(x, y) ^ (1 ^ S{x - ReE,) S{y - ImE,)) , (7) 

i 

is then obtained (upon recalling the identity ^ ^ = 7t6{x) S{y)) as 

pi^^y) = (8) 
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with z = X + iy. Here we define 



G{z)^{^irG{z)) (9) 

and < . . . > denotes, as usual, averaging with respect to the probabihty ensemble 
from which H is drawn. A small subtlety is that while G{z) is ostensibly a function 
of z only, {G{z)) depends on both z and z* . A careful discussion is given in [0]. It is 
also useful to define the Green's function Gq{z) = j^jj^ in the absence of impurities 



and the corresponding Gq{z) = jj ti ■ For the spectrum in (11), we can determine 



Gq{z) explicitly 

In the literature [0], the potential 

"^^^ ^ 2^ ^ ~ ^"^^ " ^^^^ 

has also been introduced. We see that the density of eigenvalues is given by 

which can be interpreted as the charge density generating the two-dimensional elec- 
trostatics potential ^{z,z*). See [0] for a careful discussion. 

The following probability distributions were studied in FZ and BZ: the sign dis- 
tribution 



P(^w) = ^[6{w -r) + 6{w + r)] (12) 
with some scale r, the box distribution 

P{w) = ^e{V'-w') (13) 

obtained by "smearing" the sign distribution, and the Cauchy distribution 



n«^) = -^^, (14) 



with its long tails extending to infinity. 
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We might also wish to consider the effect of diluting the impurities by setting 
randomly some fraction d of the w^'s to zero. In other words, given P{w) we can 
consider 



Pd{w) = d 6{w) + (l-d) P{w) (15) 

with < d < 1. As d increases, there should be more extended states. 

In BZ, a general formalism for studying this problem was established. A formula 
for the Green's function 



G 



1 -1 



Go 



, 1 + VkG'o 



n -1 



was obtained. Here we have defined the projection operator 

Pk = \k){k\ 



Go (16) 



(17) 



onto the site k and the repeated scattering amplitude on the impurity potential at 
site k 



Vk = 



Wk 



1 - Wk{k \Go\k) 

The reader is invited to average the trace of (|T6|) with his or her favorite P [w), and 
thence to obtain the density of eigenvalues p{x,y) by (^. For an arbitrary P{w) it is 
presumably not possible to average (|T^) explicitly and obtain {G{z)) in closed form, 
not any more than it is possible to obtain {G{z)) in closed form for the hermitean 
problem. Indeed, (^) was derived by only assuming translation invariance for Hq, 
whether hermitean or non-hermitean. It is, however, always possible to expand (p!6| ) 
to any desired powers of v and average. Some results are given in BZ. 

Remarkably, explicit results can be given for the Cauchy distribution. These 
results were also obtained independently by Goldsheid and Khoruzhenko |l^. Here 
we state the results from BZ. 

The density of eigenvalues consists of two arcs and two wings. Away from the real 
axis we have 
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p{x,y) = po{x,y + -f) 9{y) + po{x,y - -f) 9{-y) (19) 



with pq the eigenvalue density corresponding to Gq{z). In other words, the two arcs 
of the eUipse are pushed towards the real axis by a distance 7. On the two wings, the 
eigenvalue density is 



Pu,ing{x,y) = 5{y) 9{x -x„i„(7))-5^ — (20) 



where 



5(x,7,t) = V(x2+72)2 + 2t2(72_a;2) + ^4 (21) 



and 



Xminil) = {\j{t sinh hf - 72) tanh /i (22) 
The critical value of 7 at which the arcs disappear is given by 

7c = it sinh h (23) 

All states are now localized. 

In this model, 7ci is 0"*". The wings appear as soon as 7 > 0. We see that for 
insufficient non-hermit icity, h < he where 



K = sinh-i (|) , (24) 



all states are localized. Note that in accordance with Anderson et al |jT2|], for any 
finite non-zero 7, no matter how small, there are localized states for h small enough. 
Trivially, we can express (^3]) in yet another way, by saying that for a given non- 
hermiticity h and site randomness 7, we need the hopping to exceed a critical strength 



(25) 



before states become delocalized. 

For P (w) other than Cauchy, we are not able to obtain p{x, y) exphcitly. 

In FZ, it was pointed out the the problem (^) simplifies in the maximally non- 
hermitean or "one way" limit, in which the parameters in (j^) are allowed to tend to 
the (maximally non-hermitean) limit /i — > oo and t — such that 

(26) 

The particle in (^) can only hop one way. The spectrum @ changes into 

'2,7TtTi 

En = exp^, (n = 0,l,---,iV-l), (27) 

and the ellipse associated with (^ expands into the unit circle. In BZ, bounds on 
the domain over which the density of eigenvalues is non-zero were obtained for ar- 
bitrary P(w)'s with bounded support. Furthermore, it was shown that for the sign 
distribution the original unit circle spectrum of the deterministic is distorted by 
randomness into the curve 



= r' + e'\ 0<e <2tt (28) 

in the complex z plane. (See figure 2.) Clearly, = 1 is a critical value of r. For r < 1, 
the curve ( pHj) is connected, enclosing a dumb-bell region free of energy eigenvalues. 
At r = Tc, the curve traces out a figure eight. For r > 1 it breaks into two disjoint 
symmetric lobes that are located to the right and to the left of the imaginary axis. 

Incidentally, it is instructive to understand this phase transition in the language 
of (|I|) and in particular the appearance of a zero eigenvalue at r = 1. In the one way 
limit we have for r = 1 

= + w,P, (29) 

A static solution ^ = is obtained with Pj^i = — WjPj. Since Wj = ±1 occur 
with equal probability, we have a 50% chance of matching the boundary condition 
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Pn+1 = Pi- Thus, we have a zero eigenvalue in the averaged density of eigenvalues 
and hence a static state. 

Here we should mention that, strictly speaking, it is not clear that the simple 
argument mentioned above showing that localized states have real eigenvalues and 
that complex eigenvalues imply extended states would still hold in the one way limit. 
The spectrum should certainly behave smoothly since after all we are simply letting 
the clockwise hopping amplitude go to 1 and the anti-clockwise hopping amplitude 
go to 0. We have checked this numerically and can easily check this in the explicitly 
solvable Cauchy model. We have studied in FZ and BZ the localization properties 
of various models in the simplifying limit of having only one impurity present, and 
showed analytically that the states on the wings are localized. Indeed, we found that 
the localization length L{E) diverges as 

as E varies along the wings and approaches the points ±E^, at which the wings are 
attached to the arcs of complex eigenvalues. We conjectured in BZ that the behavior 
in ( PPD is generic. We have verified numerically that in the one way limit the states on 
the wings are indeed localized. We will thus take as our working assumption that the 
argument mentioned above continues to hold in the one way limit, and hence refer to 
E^: as the energy of the localization-delocalization transition in the one way limit as 
well as in the more general two way case. In the next section, we give a simple proof 
that this assumption is indeed correct. 



10 



2 Spectral Curves of Non-Hermitean Hamiltonians 



In the explicit solutions mentioned here, the eigenvalue density p{x, y) traces out a 
one-dimensional curve. This feature also appears to hold in various numerical studies 
of p{x,y). Recently, Goldsheid and Khoruzhenko |jl6|| (to be referred to as GK) gave 



a rather nice proof of this fact. We will describe a simplified version of their proof 
in the appendix. (We mention in passing that the spectra of some of the random 
hopping models studied in FZ are explicitly (and interestingly) two dimensional, but 
it can be readily checked that these models violate the assumptions made in GK.) 

We now give, in the maximally non-hermitean limit, a simple "one-line" proof that 
the spectrum of H is one dimensional. In FZ it was found by direct evaluation that in 
the maximally non-hermitean limit the determinant of {z — H) simplifies drastically 
to (we take N to be even for definiteness) 



' N \ 

det {z-H)= \ X{{z - Wi)j - 1 (31) 

(For arbitrary t and h the corresponding formula is considerably more complicated.) 
Note that (^Tj) is completely symmetric in the {wj}, and thus, for a given set of 
site energies it is independent of the way the impurities are arranged along the ring. 
(Incidentally, this property is destroyed if we include a next nearest neighbor hopping 
= s in addition to iYj^i+i = 1, and thus is not solely a consequence of the 
particle hopping one way.) 

For an alternative derivation of ( |3lD let us write down Schrodinger's equation in 
the maximally non-hermitean limit 

{E-wj)^j = 7^,-+! (32) 

N 

Starting with ipi = i/^n+i, we have ipN = {E — wn)~^ ipi-, ipN-i = D {E — 

j=N~l 

N 

Wj)~^ ijji, and so on, until finally ipi = Jl {E — Wj)^^ ipi. Hence the eigenvalue 
equation is indeed obtained by setting the right hand side of (|3TD to zero. 



With this simple form for det{z — H) we obtain the potential defined in (|10D as 
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= ^ ( I (nf=i - ^.)) - 1 T) (33) 

Let us understand the behavior of this function. Simply stated, the essence of our 
proof is to note that the quantity Hili — Wi) in ( P5D goes suddenly from being 
neglibible compared to the 1 in (^) to being completely dominant over the 1. The 
suddenness of the transition is due to the large N limit. The quantity in question is 
the product of factors, with tending to infinity. Thus, if each of the factors is 
larger than 1 on average this quantity is exponentially large, and if each of the factors 
is smaller than 1 on average, exponentially small. 

We can flesh out this argument somewhat. Consider the absolute magnitude of 
H'iLi {z — Wi) and ask if it is larger or smaller than 1 in (|5BD. For convenience, define 



Qiz) = -Vf ^og \l[ {z- Wi) I 

i=i 

1 ^ 

= JTt ^(^9 \ z - Wi \ (34) 
1=1 

Rather than dress up our proof with mathematical abstractions, we will describe 
what actually happens in specific numerical terms. To focus our thoughts, let us 
imagine we are doing numerical computation with = 10^ with P (w) the box 
distribution with V = 1, that is, the Wj's are evenly distributed between —1 and +1. 
Consider g(1.5), say, the average of 10^ numbers ranging between log 0.5 and log 2.5 
and equal to log 1.5 on the average. Obviously, q{1.5) is most likely positive. In 
contrast, consider g(0.5), the average of 10'^ numbers ranging between log \ —0.5 | 
and log 1.5. Evidently, the function q{x) is most likely negative for x less than some 
critical value x* and then most likely goes positive for x greater than x*. (For a given 
realization {wi}, the function Q{z) = e^ii^) is the absolute value of a perfectly 
smooth polynomial with zeroes at the Wj's, and does not fluctuate significantly 
from realization to realization since for A^ large, realizations in which most of the WiS 
are almost equal to +1, say, are exponentially unlikely.) In the N —>■ oo limit, the 
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function Q{x) jumps from zero for x < to infinity for x > x*. To see that 
this jump is abrupt in the large N hmit, consider that in going from = lOHo 

= 10^, we basically take the Q{x) for = 10^ and raise it to the power 10^. 
Similarly, for a fixed 6, as we increase r there is a critical r^,{6) at which Q{re^^) jumps 
from zero to infinity. 

From (^) we see that for a given and for a given realization, ^{z,z*) goes 
from essentially zero (when Q{z) = \ YiiLi {z — Wi) \ is small compared to 1) to 
a non-zero value (when Q{z) is large and the 1 inside the logarithm is (^) can be 
neglected.) The transition occurs when Q{z) = 1 and hence q{z) = 0. As A^ ^ oo, 
the transition region becomes narrower and narrower and so the transition region is 
one-dimensional . 

Hence we see that in the complex plane there is a curve across which ^{z,z*) 
changes from 

^{z,z*) = (35) 

to 

Hz,z*) = {^log\Ul,{^-w^)\') 

= { log I z — w I) (36) 

We have proved that the spectrum in the maximally non-hermitean limit is one- 



dimensional and that the spectral curve is given by the dividing line between (|35| ) 
and (M 



{log \z-w\) = (37) 

Our proof can be made more mathematically rigorous by the usual considerations 
that for A^ large the properties of H for a given realization {wi} are the same as the 
averaged properties of H for N = oo. 
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To summarize, we want to emphasize how simple our proof is. First, we note that 
in the maximally non-hermitean limit the determinant of {z — H) is given simply by 
(Pip. Then we just look at (^Tj) and ask which of the two terms on its right hand side 
is larger. The transition between the two regimes gives the spectral curve. 

We are now in a position to understand the wings in an extraordinarily simple 
fashion. In the region where (|36|) holds the density of eigenvalues is given, according 
to ([Tl|) by simply 



Pwing {x,y) = 6{y) 6{x ~ x^)P{x) {31 



Here x^, is defined by the region where (|36D holds. We see that the wings stick out 
from the central "elliptical" region, since inside the curve defined by (^) the potential 
$ vanishes and hence p vanishes. 

This discussion also shows that the density of eigenvalues on the wings simply fol- 
lows the probability distribution of the site energies. We can check this result against 
our explicit result in (pO|) . In the one way limit, t — > and hence B x^ + 7^- 
We see that ( pOD indeed reduces to (pS]). We also understand why there are no wings 
in the sign model in the one way limit. 

While (38) is not immediately obvious from the Schrodinger equation (32), we 
can see heuristically from (^1]) that for a given W}^ how one might be able to find an 
eigenvalue near by writing E = ^ for z in (^T]). The condition for a 
solution is that the product V^i^^kh^k — Wi) is approximately equal to ^. According 
to (|38D , this is possible in a statistical sense for Wk large enough (larger than x^,.) We 
see from (R^) that the solution is clearly peaked on the k^^ site. 



The reader who has read FZ and BZ may recognize that the way we derived (|37|) 
was essentially how the spectral curve (^) for the sign model was obtained. In the 
sign model, each Wi is equal to either +r or — r with equal probability. For large, 
in a given realization, most likely about half of the Wj's are equal to +r, and the 
rest equal to — r. Thus Q{z) = \ T^fLl {z — Wi) \ ~ \ z'^ — r'^ \^ . Thus, for 
> 1 the 1 in ( p3D is completely overwhelmed, while for \ z^ — \ < 1, 



^2 _ ^2 
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the 1 dominates. The spectral curve is given by | 2;^ — | = 1 as in 

We also mention that the bounds in BZ, alluded to in section I, were obtained by 

a related argument based on whether | {(^jz^^ ) I is less than 1 for some integer /. 
The proof of Goldsheid and Khoruzhenko is more involved but also more general, 

and applies to the non- maximally non-hermitean limit. In the appendix we give a 

much simplified version of their proof, purging what we do not need for our purposes. 
Finally, the discussion in this section enables us to calculate the localization length 

L {E) on the wings. According to the discussion following (|32D, if we move n lattice 

spacings away from the site where the wave function attains its peak value, the 

absolute value of the wave function | ip \ decreases by a factor 



n - Wj 



e 



-E log \E-Wj\ 



g-n (i E log \E-Wj\) ^gg^ 



where the product and sum in (^) involve n terms. For n large, we see that the 
localization length L{E) is none other than the inverse of the function q{E) defined 
in (H): 

^i^)-^) (40) 
Indeed, q{E) is positive on the wings. As E approaches E^ (also called in some 
contexts) at which the wings are attached to the arcs of complex eigenvalues, q{E) 
vanishes according to (|35D and hence generically 

HE) ~ ^^4^ (41) 

as we conjectured. 
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3 Critical Points in Maximally Non-Hermitean Mod- 
els 



We will try to extract as much as we can from {p7\) for a general probability distri- 
bution P{w) of the site energies. (For simplicity, we will, as before, take P{w) to be 
even.) Later, we will obtain explicit results for specific choices of P{w). The equation 
flSTl) which we write as 



( log ((x - wf + y^)) = Q (42) 

defines a curve s{x, y) = in the complex plane. With P {w) even, we have 
s{x, y) = s{—x, y) and so we will henceforth focus on the first quadrant x > , y > 0. 

Let us focus on the point y^: at which the spectral curve s intersects the imaginary 
axis. In general, the probability distribution P (w) is characterized by one or more 
parameters 7 which measure the randomness. We will by convention take 7 = 
to describe no randomness, and associate increasing randomness with increasing 7. 
Define 

/(l/,7) = {log iy' + w') ) (43) 
The desired function ?/*(7) is determined by 

/(2/*(7), 7) = (44) 
We know that with no randomness y*(0) = 1, and we expect that as 7 increases, 
?/*(7) decreases until a critical point 7c at which point ?/*(7c) = 0. 
These expectations are easy to prove. First, note 

I- = 2l/(^^) > (45) 
ay y^ + 

(recall that we restrict ourselves to the first quadrant) and hence fiy,^) is a mono- 
tonically increasing function of y. At y = 0, the function starts with the value 
/(0, 7) = ( log w'^), and increases to / ~ log y"^ for large y. Thus, if ( log w"^ ) < 0, 
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the function f{y,'y) vanishes at some y*{j). But if ( log w"^ ) > 0, there is no y*{j). 
Hence, the critical value 7c is determined simply by 

{log w"^) = (46) 

A critical transition occurs when the average of the logarithm of the site energy 
squared vanishes. 

Let us give a specific example, the box distribution in (|13[). Note for our purposes 
an overall multiplicative factor in /(?/, 7) is irrelevant and thus we can just as well 
write 



rV 

f{y,l) = / dwlog{y^ + w^) 
Jo 

= 2y arc tan ^ + ^ (-2 + log{y^ + V^)) (47) 

(Here the generic parameter 7 is denoted by V.) In particular, up to an irrelevant 
overall factor, we may write (^6]) as 

-2 + log = (48) 

and thus we find 

Vc = e (49) 

As we increase the randomness (as measure by 7), we expect the eigenvalues to 
migrate towards the real axis, and hence y^Xl) to decrease. This expectation can be 
easily proved for a general P (w). To avoid notational clutter, let P (w) contain only 
one parameter 7 relative to which the random site energy w is measured so that we 
can write P (w) = 0(7) p(^). (Thus, for example, for the box distribution ([TsD , 
the one parameter 7 is called V, and 0(7) = and p(^) = 6 (l — (■^)^).) 
Differentiating the defining equation for ?/*(7) 
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roc yj 

/ dw Pi-) log {yl + w^) = Q (50) 
JO 7 



7 

with respect to 7 and integrating by parts, we obtain 



thus proving that ^ < 0, as desired. 



Indeed, it is easy to obtain a general result for small randomness. Expanding (H3) 



{log (yl + w^)) = Q (52) 
for y^. = \ — 6y, we obtain {—26y + w'^) = and hence immediately 

y*ii) = 1 - l{w') + --- (53) 

for small randomness. In particular, for the box distribution, y^iV) = 1 — | + ■ ■ ■. 
Also, for the box distribution, in the other extreme, when randomness reaches the 
critical value Vc given in (^), y*{V) ~ ^i^c ~ V) vanishes linearly. 

It is now just a matter of doing the appropriate integral to obtain 7^ for various 
probability distributions of the reader's choice. Notice that in many cases we may be 
able to determine 7c analytically without being able to determine 1/(7) analytically. 
Let us give some examples. The simplest is the sign distribution (0) for which 
fiVy ^) = ^og (y^ + r^) and thus 



y^ = Vl-r2 (54) 

giving Tc = 1 in agreement with (^). It is also interesting to look at the Cauchy 
distribution (W^ for which we find 



°° log w"^ TT 

^2 + 7^ 7 



dw z::^ — 2 = tJ^9 7 (55) 



and thus 
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7c = 1 (56) 

in agreement with (p3|). 

A class of probability distributions for which {log w"^) is easily evaluated is the 
generalization of the box distribution defined by P (w) oc \ w \^ 6 (V"^ — w"^). We 
obtain 

Vc = (57) 

and thus recover (|49|) for k = 0. Numerically, we have = (2.718, 1.649, 1.395, ■ ■ ■) 
for k = (0, 1, 2, ■ ■ ■). The general trend agrees with our intuition: for a given V, as k 
increases, the random site energies become larger, and hence Vc should decrease. As 
— ^ oo we have = 1, in agreement with = 1 for the sign model as expected. 
Similarly, we can study the "parabolic" distribution P {w) oc (V^ — w"^) OiV"^ — w"^) 
and distributions similar to it. We find 

Vc = = 3.794 (58) 

-^^ 

For the Gaussian distribution P (w) oc e we obtain 

= 2e^ = 2.669 (59) 

where denotes Euler's constant 0.577. For the exponential distribution P (w) oc 
we obtain 

ttc = = 1.781 (60) 

For the same a, the exponential distribution is "more random" than the Gaussian 
distribution, and hence we expect Oc to be smaller for the exponential than for the 
Gaussian distribution. 

We have obtained the value of 7c, the critical measure of randomness, for a variety 
of probability distributions P{w). What happens when the randomness reaches its 
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critical strength? In the Cauchy model, we know that at 7c all eigenvalues become 
real and the extended states disappear. In general, with increasing randomness, 
eigenvalues tend to drift towards the real axis, and thus we might guess that at 7c all 
eigenvalues become real, as in the Cauchy model. This however is not true in general. 
In fact, we have already encountered a counter-example, the sign model. According 
to (|28D, for r > Tc = 1, the spectrum splits into two disjoint lobes. It is indeed true, 
however, that y*(r) goes to zero as r increases towards Vc, as indicated by (^4|). 

This instructive example tells us that to determine whether the eigenvalues have 
all collapsed onto the real axis, it is not enough to study the intersection of the 
spectral curve with the imaginary axis. We have to study the intersection of the 
spectral curve with the real axis. 

According to (^2l), the point Xj.l'-f) at which the spectral curve intersects the real 
axis is determined by 

( log ix,-wf) = (61) 
With P {w) even, we can symmetrize in w and define 

g{x,'y) = {logix^-w^) (62) 

From (^21) directly, or from (^), we see that g{x,'y) is given by the analytic contin- 
uation of 2/(?/,7). The sign flip between f{y,'y) and g{x,'y) makes a difference of 
course: it is no longer true that g{x,'~f) is a monotonically increasing function of x. 

However, we will now prove that if P (w) does not increase for w positive and 
increasing, then g{x,'^) is indeed monotonically increasing as x increases. The proof 
is straightforward. Using P (w) > P{x) ioi w < x and —P {w) > — P{x) for 
w > X, we find that 

' > Pix)r dw^l-^ = (63) 



Ax dx Jo x"^ — w'^ 

(where the last integral is understood in the principal value sense, of course.) For 
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large x, g{x,'y) ~ log x'^. Thus, if 5'(0,7) = {log w'^) is negative there exists a 
solution X* (7) to the equation 

( log {xl - w' f) = (64) 

When 5f(0,7) goes positive, there is no solution, and thus the critical value of ran- 
domness is given by 

( log w^) = (65) 

in agreement with (^6|), as it should. 

The behavior of a;* (7) for small randomness is also easy to extract. Writing 
x^i^j) = 1 + 6x and expanding (|6^), we find immediately that 

xM = ^ + I + (66) 

(This result and ( [5^ ) are of course in agreement with the appendix in FZ.) That x*(7) 
vanishes as 7 imply that a;* (7) has a maximum at some critical value 7*. 

1 . 4 
1.2 
1 

. 8 
. 6 
. 4 
0.2 

0.5 1 1.5 2 2.5 

Figure 3: The energy at which the localization-delocalization transition occurs, a;*(\^), 
as a funtion of V for the box distribution. It has a maximum at K ~ 1.81 and 
vanishes at Vc = e. 
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We show the behavior of x^{'-f) as 7 varies for the box distribution in fig(3). We 
can evaluate g{x, V) exactly to find (up to irrelevant overall factors) 

g{x,V) = {V - x) log {V ~ xf + {V + x) log {V + xf - 4V (67) 

(where as before we remind the reader that the generic parameter 7 is called V for the 
box distribution.) To determine 7* (that is, K in this context), we simply differentiate 
g{x, V) to obtain 

- <{V.) = 1 (68) 

Inserting this into g{x^:(y*), V*) = we find readily that K = cosh 9^ ~ 1.81 and 
x^,(K) — 1-51 where 6^, is the solution of the equation 6^: tank 9^ = 1. 

For V Vc = e we find x^,{V) ^ y^2Vc (K — V). We can evidently define a 
variety of critical exponents here: for example, 

X, {V) ~ (K - Vy (69) 

The alert reader may have noticed that for the Cauchy distribution x*(7) = a/1 — 7^ 
and instead of increasing with 7 as shown in ( |66D x*(7) decreases with 7. This is in 
fact consistent with (|66|) since the second moment {w"^) does not exist for the Cauchy 
distribution. The critical exponent P again takes the generic value /? = |. 

For a P (w) which increases as w positive increases, the behavior we just sketched 
for x=k(7) fails. Our counter-example is the sign model, in which 

g{x,r) = log {x^ - r'^f (70) 

a function which starts at log r^, goes to —00 to x = r, and then rises from —00. 
Thus, at Tc = 1 we go from a single solution for x*(7) to two. (As said before, we 
always work in the first quadrant.) This is consistent with (pSl). 

In summary, for P {w) which does not increase as w positive increases and which 
has a second moment, x,,(7) increases as the randomness 7 increases, reach a max- 
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imum at 7*, and then decreases towards zero as 7 — > 7c. Meanwhile, ?/*(7) de- 
creases monotonically towards zero as 7 increases towards 7c. Note that x^{'y) has 
the important physical meaning as the critical energy E^: (7) at which the localization- 
delocalization transition occurs for a given measure of randomness 7. 

Let us use the terminology the "arc part" to refer to the spectrum of H consisting 
of complex eigenvalues. Thus, increasing randomness squashes the arc part towards 
the real axis, elongating it in the x direction, until the critical randomness 7=,,, after 
which it shrinks steadily, giving up spectral weight to the wings, until it disappears 
altogether at 7c, at which point all states become localized. 

Using the sort of techniques developed here, the reader can readily deduce further 
properties of the spectral curve. For instance, with P (w) even and writing z = re*^, 
one has 

{log \ {z-w){z + w) \ ) = {log {r"^ - 2 r"^ w"^ cos 29 + w^) ) = (71) 

For 6* > ^ in the first quadrant, cos 29 is negative, and one can deduce that ^ < 
upon differentiating (ffll). 
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4 Appearance and Disappearance of Wings 



Starting with the spectrum with no disorder present, we expect that as we increase 
the randomness wings appear at some critical value of the randomness 7ci, as defined 
earlier. For the Cauchy distribution we saw explicitly that 7ci = 0+; in other words, 
wings appear as soon as the random site energies are switched on. As we will see, 
this is due to the fact that the Cauchy distribution has infinitely long tails. 

In this section we consider P (w) with support up to Wmax = 7- (For example, 
for the box distribution the maximum value of w is V, which is what we denote 
the generic parameter 7 by for this distribution.) According to (0), the right wing 
extends out to Wmai, = 7. On the other hand, the right wing starts at x^, (7). Clearly, 
the wings disappear when ^^(7) becomes larger than 7. Thus, we conclude that for 
this class of distributions •jci is determined by 



Again assuming as before that P {w) contains only one parameter 7, we see that 7ci 
is the solution of 




(72) 




(73) 



Changing variable in ([73| ) by w = 7M, we obtain the result 






(74) 



which is to be compared to 



Jq du p{u) log 




(75) 



/o dupiu) 



In particular, for the box distribution we have 



1 



(76) 



= -e 
2 
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compared to 



(77) 



For the "parabolic" distribution mentioned earlier we have 



1 5 



(78) 




compared to 



(79) 



Notice that the argument shows that if tiimax = cxd, as is the case for the Cauchy 
distribution, then the right wing extends out to infinity, and the wings are always 
present as long as 7 > 0. Thus, 7ci = 0"*" as we have already learned for the Cauchy 
distribution. Physically, for P {w) with Wmax — 00 the site energies have a non-zero 
probability of attaining arbitrarily large values and a small amount of randomness 
can cause the wings to appear. 

It is easy to show that if P (w) does not increase for w positive and increasing, 
then 7^^ < 7^. Start by noting that (1 — < (1 — m^) for < m < 1, and 
that for the class of P(w)'s under consideration (— log v^) > (— log (1 — u)^)- 
Combining these two inequalities we find the desired result. 
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5 Spectral Curve for the Box Distribution 



In the preceding section, we studied the points x*(7) and ?/*(7) at which the spectral 
curve s{x,y) = intersects the real and imaginary axes for general distributions 
P (w) of the site energies. In this section, we determine the entire spectral curve 
s{x,y) = for the box distribution ([T3|) . 
First, we evaluate fBB) to obtain 



1 / , (m+ 

2 1 7^73iy«-i 



Hz = Vu,V) = - { log I ) -{—^ V'\\ - 1 (80) 



where we scale z = Vu for convenience. Thus, the spectral curve is determined by 
the interesting equation 



iu + 1)"+^ , 



Let us write u = s + it and define the left hand side of (|ST|) as b{s,t). We find 



^1) 



b{s,t) = {{s + lf +t^y'' {{s-lf +t^y~'' e-'""'"'''^^ (82) 
The spectral curve (with z = x + iy ) is then defined by 

In particular, x*(7), the critical energy at which localization-delocalization tran- 
sition occurs, is determined by the remarkable function 

6(.,0) = (s'-l) C-±^r (84) 
s — 1 

In solving (|83|) it is essential to pick the correct branch of the arc tan function, which 
is defined up to multiples of vr. Let us take the conventional definition of arc tan as 
a function ranging between — | and |. Then we see that for < i _ b{s,t) has 
to be multiplied by a factor of e'^*. This ensures that b{s,t) is continuous as t passes 



through a/I — s^. We show in figure (4) the solution of ( ^31) for V = 1. The spectral 
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Figure 4: Spectral Curve for the box distribution with V = 1, plotted in the first 
quadrant 

curve is "elliptical" in shape. We find numerically that as V approaches Vc = e the 
fluctuation between runs even for = 200 to be quite significant. These fluctuations 
are described by the density-denstity correlation function discussed in BZ. 

We showed in figure (3) the behavior of the energy of the localization-delocalization 
transition defined by 

biy^^) = ^ (85) 

The function behaves as discussed in the last section. In particular, at = 1 

the solution of 

ix' - 1) = (86) 

X — 1 

is given by a; = 1.17697. Note also that 6(1,0) = 4 and so for V = |, we have 
x^(y) = V. Recall from the preceding section that Vd for this model is precisely |. 
We see that indeed the wings disappear at this critical value of the randomness. 

Given the complexity of (|83|) it would be a challenge to determine the spectral 
curve in the non- maximally non-hermitean case. 
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6 Dilution 



As mentioned in the Introduction, the effects of dilution were considered in BZ. We 
point out here that the effects of dilution can be quite drastic. The sign model is 
most amenable to analytic study with (p!5|) 

Pd{w) = d5{w) + \{^-d) {5{w-r) + 5{w + r)) (87) 
The spectral curve is determined by 



2\d 



\x — y — r \ + 4:X y 



Consider for simplicity the equation for y^{r). For no dilution (i = 0, this becomes 

yl + = 1 (89) 
as in (|69|). For any amount of dilution, c? > 0, we have 

yl\yl + r'Y"' = 1 (90) 



Note the left hand side of both (|89|) and (pO]) are monotonic in y^,, but the left hand 



side of (pO|) starts at instead of and so a solution of (|90|) for y^, always exists for 
any > 0. This is in sharp contrast with (|89|), for which y^: exists only for r < 1. 
As shown in fig(5), for r > 1 and any (i > a new "lobe" centered at the origin 
appears in the spectrum, with this lobe intersecting the imaginary axis at y^{r,d). 
We can easily solve ( pOD for c/ = |, obtaining 



1. 1 



y.{r,d = -) = W-(vW4 - r2) (91) 



More generally, for any rf, we see that for small r 

1 

y^[r,d) 

and for large r 



y,{r,d) ^ - (l-c/)r2 + ... (92) 
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Figure 5: The diluted sign model with d = 0.1 and r = 1.2 

We see that there is a critical transition with dc = 0, as emphasized in particular 
by®. 

More generally, recalling the definition of f{y,'y) = { log {y^ + w"^)) from (|43|) 
associated with P(w), we have associated with Pd(w) the function 



f{yn,d) = dlogy' + {1 - d) fiy,j) (94) 

with the intersection of the spectral curve with the imaginary axis y^,{'~f,d) given by 
/(y*, 7, (i) = as in (|4^). Evidently, as the dilution goes to its maximum value — > 1, 
y* 1- 

For the Cauchy distribution, we have 



fiy,l) = log {y + if 
With no dilution, ?/*(7) = 1 — 7 as in FZ. For small rf, 
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(95) 



y,{^,d) = 1-7 - {log{l-^)) d + ■■■ (96) 
Clearly as 7 nears its critical value 1, a small dilution has a big effect. For c? ~ 1, 



y,(7,d) = 1 - ( W(l + 7)) + ••• (97) 

Further study of the effects of dilution would be interesting. 
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7 Appendix: The Proof of Goldsheid and 
Khoruzhenko Simplified 



Here we give a simplified version of tlie proof of Goldshied and Khoruzhenko. Let 
us go back to (|^). The Hamiltonian H would be hermitean had the entries Hj^i and 
i^iAT been absent. It is natural to define the hermitean Hamiltonian by 

H = Hh + V (98) 
where all entries of V are zero except 

V,N = le-'^' (99) 
(At first sight, the transformation of H into H appears to be a potentially dangerous 



manuever [|T7| as H is an N hy N matix with an entry which explodes exponentially 
as ^ oo; in fact, we are going to exploit precisely this feature.) 

Since H and H have the same spectrum, we can obtain the desired density of 



eigenvalues of H by 



with 



^{z,z*) = ^Alog \ det {z - H) \^) (101) 



2N 



In analogy to (|101|) we also define 



^h{z,z*) = ^{log \ det {z - Hh) \^) (102) 



2N 

corresponding to the hermitean Hamiltonian Hh. Using 



[z-H) = {I-VG) {z-Hh) (103) 
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with 



G . ^ (104) 

(not to be confused with tlic G we defined in (|]); we use the same symbol to avoid 
notational clutter) we have immediately 

^z,z*) = Mz,zl + \det{I-VG)\^) (105) 

Since V has only two non-zero entries, the determinant in (|105|) can be evaluated 
immediately to be 



det (/ -VG) = {1 - ^- e-''^GNi) (1 - ^e^'^d^) - ^ Gu Gnn (106) 

The large behavior of det {I — VG) thus depends on whether Gin vanishes faster 
than e~^^ or not. 

If Gin vanishes faster than e~^^, det{I — VG) is bounded. (Since ;^trG is the 
Green's function of a hermitean Hamiltonian it has a finite large limit, and Gu 
and Gnn are expected to be the same order of magnitude as ;^trG.) Then from ( |105|) 
we have 

^z,z*) = Mz,z*) (107) 
In contrast, if Gtm does not vanish faster than e~^'^, then 



det {I - VG) - ^ e^'' Gin (108) 



as ^ oo, and we have to keep the second term in ( |105| ). From the definition (^ 
and the form of we obtain immediately 



and hence from (|102|) 
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{log I Gin \) 



log 
log 



Inserting (|11CI|) in ( |108|) we obtain 



2) 

2. 



{log I det {z — Hh) \) 



(110) 



l-{log\det{I-VG)\) = 1 



log ( ^e' 



- Ar$;,(z,z 



<l>h{z,z*) 



;iii) 



At this point we might as well set t = 2. Thus, combining (|105|) and ( |111| ) we see 
that if Gin does not vanish too fast we have 



^z,z*) = h (112) 
According to ( |100| ), H and hence H do not have any eigenvalues in this region. 



We thus obtain the result of GK that the dividing line between (|107|) and (Hi 
is the curve defined by 



h 



(113) 



This forms the "arcs" of the spectral curve of H. The wings are determined by the 
solution of ( |1071 ) which lies outside the curve defined by (|113|) . 

The defining equation ( |113| ) for the spectral curve, while elegantly compact, can- 
not be explicitly solved in general. We need an exact analytic expression for $/i(z, z*) 
which amounts to an exact determination of the density of eigenvalues of the her- 
mitean Hamiltonian Hh- To the best of our knowledge, this is available only for the 



Cauchy distribution ([T^) . Indeed, Goldsheid and Khoruzhenko [|T^ are then led to the 
curve in (|T^, which was derived in BZ by a different method. 

Thus, in order to obtain explicit expressions for the spectral curves for various 
P(w)'s, we are compelled to got to the maximally non-hermitean or one way limit 



33 



(|26|). Going back to ( |1 1 1| ) we see that in this hmit the spectral curve is actually 



defined by setting h = in (|113|) . 

Next we observe another drastic simplification. In this limit, the hermitean Hamil- 
tonian Hh goes over to a purely diagonal matrix with the j*'* diagonal element equal 
to Wj, and thus the function ^h{z,z*) defined in ( |102| ) simplifies drastically to 

Mz,z*) = ^ {log \z-w\^) (114) 

Thus, in the maximally non- hermitean limit, the spectral curve of the non- hermitean 
Hamiltonian H in is determined by the remarkably compact equation 

{log \ z-w\^) = (115) 



in complete agreement with (^) and our simple "one-line" proof. It is interesting to 
note that in this respect the non-hermitean problem is again simpler than the her- 
mitean problem. There is no analogous simplifying limit we can take in the hermitean 
localization problem. 



34 



References 



[1] F. Haake, F. Izrailev, N. Lehmann, D. Saher and H. J. Sommers, Zeit. Phys. B 
88 (1992) 359. 

H. J. Sommers, A. Crisanti, H. Sompolinsky and Y. Stein, Phys. Rev. Lett. 60 
(1988) 1895. 

M. A. Stephanov, Phys. Rev. Lett. 76 (1996) 4472 . 



R. A. Janik, M. A. Nowak, G. Papp and I. Zahed, preprints |cond-mat/9612240 
|iep-pli/960632-^ 



R. A. Janik, M. A. Nowak, G. Papp, J. Wambach and I. Zahed, preprint |hep- 



ph/9609491 



E. Gudowska-Nowak, G. Papp and J. Brickmann, preprint |cond-mat / 9 70 1 1 87 



Y. V. Fyodorov, B. A. Khoruzhenko and H. J. Sommers, Phys. Lett. A 226 
(1997) 46; Phys. Rev. Lett. 79 (1997) 557; B. A. Khoruzhenko Jour. Phys. A. 
(Math. Gen.) 29 (1996) L165. 

M. A. Halasz, A. D. Jackson and J. J. M. Verbaarschot, preprints, jiep- 
lat/970300q , |hep-lat/9704007 . 



J. Feinberg and A. Zee, preprint |cond- mat / 9 704 1 9l| , Nucl. Phys. B 501 (1997) 
643. 

[2] J. Miller and Z. Jane Wang, Phys. Rev. Lett. 76 (1996) 1461 
J. T. Chalker and Z. Jane Wang, preprint |cond- mat / 9 704 1 98 . 



[3] D.R. Nelson and N.M. Shnerb, |cond-mat/9708071 



[4] N. Hatano and D. R. Nelson, Phys. Rev. Lett. 77 (1997) 570; for a more detailed 
discussion see N. Hatano and D. R. Nelson, preprint |cond-mat / 9705290 . 



[5] One of us, A. Zee, thanks P. Le Doussal for a discussion of this point. 
[6] K. B. Efetov, preprint |cond-mat/ 970209 1| , 9706055. 



35 



[7] J. Feinberg and A. Zee, preprint |cond-mat/9703087| , Nucl. Phys. B 504 (1997) 
579. 



[8] R. A. Janik, M. A. Nowak, G. Papp J. Wambach and I. Zahed, preprint |cond- 
mat/ 9705098 . 



[9] P. W. Brouwer, P. G. Silvestrov and C. W. J. Beenakker, preprint |cond- 



mat/9705186 



[10] J. Feinberg and A. Zee, preprint |cond-mat/97062T8 



[11] E Brezin and A. Zee, |cond-mat/9708029 



[12] P. W. Anderson, Phys. Rev. 109 (1958) 1492; 

R. Abou-Chacra, P. W. Anderson and D. J. Thouless, Jour. Phys. C6 (1973) 
1734; R. Abou-Chacra and D. J. Thouless, Jour. Phys. C7 (1974) 65; 
E. Abrahams, P. W. Anderson, D. C. Licciardelo and T. V. Ramakrishnan, Phys. 
Rev. Lett. 42 (1979) 673; 

E. Abrahams, P. W. Anderson, D. S. Fisher and D. J. Thouless, Phys. Rev. B 
22 (1980) 3519; P. W. Anderson, Phys. Rev. B 23 (1981) 4828. 

[13] P. Le Doussal, unpublished 

[14] A version of this argument was givenby P.G. de Gennes, J. Stat. Phys. 12 (1975) 
463, in the context of the Fokker-Planck equation. 

[15] L.W. Chen, L. Balents, M.P.A. Fisher, and M.C. Marchetti, Phys. Rev. B54 
(1996) 12798 

[16] l.Y. Goldsheid and B.A. Khoruzhenko, preprint |cond-mat/9707230| 
[17] One of us, A. Zee, thanks T. Banks for a discussion of this point. 



36 



